Is the Radio Core of BL Lac Precessing? 



R. L. Mutel 

Dept. Physics & Astronomy, University of Iowa, Iowa City lA 52242 

robert-mutelOuiowa . edu 
G. R. Denn 
Sweet Briar College, Sweet Briar, VA 24595 
gdennOsbc . edu 

ABSTRACT 

Stirling et al. (2003) reported the discovery of a 2.3 year periodic variation 
in the structural position angle of the parsec-scale radio core in the blazar BL 
Lac . We searched for independent confirmation of this periodic behavior using 
43 GHz images of the radio core during ten epochs overlapping those of Stirling 
et al. Our maps are consistent with several periodicities, including one near the 
period reported by Stirling et al. By comparing our position angle measurements 
with those of Stirling et al., we find strong consistent evidence for position angle 
variations of the inner core during the observed epochs. However, the claim 
for periodic variation is not convincing, especially when the most recent epochs 
(2000.60 - 2003.78) are included. A definitive resolution will require continued 
monitoring of the core structure over several periods. 

Subject headings: VLBI, Active galaxies, precession, relativistic jets 



1. Introduction 

There is growing evidence that at least some active galactic nuclei (AGN) display peri- 
odic behavior on a timescale of years or decades. The observational evidence consists largely 
of periodicities in optical or radio light curves or periodic variations in the parsec-scale radio 
cores mapped by VLBI. Perhaps the best documented case study is that of the BL Lac 
object OJ287 for which a period of 11.6 year was found in the visible light curve (Sillanpaa 
et al. 1988). This periodicity was confirmed by observations of a predicted outburst in 1994 
(Sillanpaa et al. 1996; Pursimo et al. 2000) but with a one year delay between the optical 
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and radio peaks (Valtonen, Lehto, & Pietila 1999). Katz (1997) suggested that the observed 
periodicity is due to the sweeping of a relativistic beam ahgned normal to a precessing ac- 
cretion disk surrounding a supermassive black hole. Valtaoja et al. (2000) suggest a binary 
black hole model in which the secondary orbit penetrates the accretion disk of the primary, 
causing a thermal (optical) pulse, followed by enhanced accretion and increased particle flux 
in the relativistic jet, which is responsible for the time-delayed radio flare. 

Precessing relativistic beams have also been invoked to explain the undulating jet struc- 
tures in several compact radio sources, including 1928-1-738 (Roos, Kaastra, & Hummel 1993; 
Murphy Preston, & Hirabayashi 2003), 0153-^744 (Hummel et al. 1997), and 4C12.50 (Lister 
et al. 2003). Britzen et al. (2000, 2001) interpreted the strongly bent relativistic jet in PKS 
0420-014 as evidence for a binary black hole system which precesses on a timescale ~ 10 yr. 
Abraham & Romero (1999) analyzed eight superluminal jet components of the parsec-scale 
radio quasar 3C273 observed over the last 30 yr. They suggest that the observed velocities 
and position angles can be interpreted by ejection of synchrotron-radiating shocks moving 
outward in ballistic trajectories from the base of a precessing inner jet whose period is 16 yr. 

These observations provide post facto evidence for precession in the central engines of 
AGN i.e., the observed jet morphology is interpreted as a record of the putative jet motion, 
but the jet oscillation itself has not been observed. However, in a recent paper Stirling et 
al. (2003) report that they have directly detected periodic changes in the parsec-scale radio 
jet orientation of the eponymous AGN BL Lacertae based on analysis of two independent 
datasets: Periodic variations in the polarization position angle at 1 mm wavelength, and in 
the direction of the innermost radio core component in 43 GHz VLBI maps. In the latter 
case, they find strong evidence of periodicity with a timescale of 2.29±0.35 yrs and an angular 
amplitude of 24.4° ± 16°. The VLBI maps spanned 23 epochs over the time range 1997.58 
- 2001.28. They modeled the source brightness at each epoch using a variable number of 
elliptical Gaussian components, including a two closely-spaced (~ 0.1 mas) circular Gaussian 
subcomponents in the core. The model component parameters were adjusted for a best-fit to 
the observed brightness distribution. They found a periodic variation in the core 'structural 
position angle' (SPA), defined as the relative position angle between the two closely spaced 
core subcomponents. 

This paper reports results from a series of 43 GHz VLBI observations of BL Lac using 
very similar observing parameters during an overlapping time range (1998.73 - 2003.82). 
Although the primary purpose of the observations was rather different (monitoring the mag- 
netic field structure of the radio jet, Mutel & Denn, 2004 in prep.), the resulting maps are 
of comparable quality to those of Stirling et al. (2003) and provide an independent test of 
the precession hypothesis. 
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2. Observations and Data Reduction 

The observations were performed at regular ~0.3 year intervals during nine epochs 
between 1998.76 and 2002.05 and a tenth epoch at 2003.82 using the 10-element Very Long 
Baseline Array ^ (VLBA). The observations consisted of a series of short scans at each of 
three observing frequencies (15.4, 22.2, 43.2 GHz) of BL Lac and three other AGN sources in 
sequence for 12 hours. We recorded both left and right circular polarization at each telescope 
and correlated all four cross-correlations pairs (RR, LL, RL, and LR) for each basehne. The 
data were recorded in standard VLBA mode using 128 MHz bandwidth and 1-bit sampling. 

Prior to self-calibration, all fringe-fitting and visibility calibration was done using the 
Astronomical Imaging Processing System (AIPS, van Moorsel, Kemball, & Greisen 1996). 
We performed amplitude calibration using the AIPS automated calibration transfer system 
(Ulvestad 1999). Fringe fitting of the parallel (RR, LL) and cross-hand (RL, LR) polarization 

correlations, as well as removal of instrumental polarization were done using the scheme 
described in (Denn, Mutel, and Marscher 2000). Hybrid mapping was done in AIPS and 
with the Caltech hybrid mapping program DIFMAP (Shepherd 1997). The resulting Stokes 
I maps had typical RMS noise level ~ 0.8 mJy per beam at 43 GHz. 



^The National Radio Astronomy Observatory is operated by the National Radio Astronomy Observatory 
by Associated Universities, Inc., under cooperative agreement with the National Science Foundation 
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Fig. 1. — 43 GHz contour maps of the BL Lac's radio core at epochs 1999.12 (left inset, 
Stirhng et al. 2003) and 1999.16 (right, present data). Contour levels increment in factors of 
two starting at 3 mJy in both maps. The overlaid crosses and circles with crosses indicates 
the size of fitted elliptical Gaussian components. 
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Figure 1 shows a comparison of 43 GHz contour maps of BL Lac observed by (Stirling 
et al. 2003) at epoch 1999.12 and our map 14 days later (epoch 1999.16). In both maps, 
eleven component best-fit elliptical Gaussian brightness models are shown overlaid by crosses 
and circles. The maps are very similar, both in dynamic range and angular resolution, as 
might be expected given the similarity in observing equipment and close proximity in time. 
However, there are small but statistically significant differences in the Gaussian component 
locations and sizes, which might reflect real differences in structure over 14 days, but could 
also be partially explained by subtle differences in the model fitting algorithms used- (Bran- 
deis VLBP package, Roberts, Gabuzda, & Wardle 1987) versus (DIFMAP, Shepherd 1997). 
Nevertheless, the close similarity in image quality and resolution of the two datasets confirms 
that the present data comprise a comparable but independent test of the claim for periodic 
behavior in the inner core. 

Gaussian component modeling of the inner core structure is illustrated in Figure 2. 
Figure 2o shows the full 43 GHz hybrid map of BL Lac at epoch 1999.73, while Figure 2b is 
a magnified view of the core only. The innermost part of the core was fit with three circular 
Gaussian components labeled CI, C2, and C3. The structural position angle (SPA) is defined 
as the position angle between the two innermost components CI and C2. It is this orientation 
angle that Stirling et al. (2003) find exhibits periodic variation which they interpret as 
evidence for core precession. Note that the C1-C2 separation is only ~ 20% of the restoring 
beam dimension along the structural position angle (lower left ellipse). Hence, the inner 
core substructure is not well resolved, which makes the Gaussian component decomposition 
non-unique. To check the reliability of the SPA determination, both Stirhng et al. (2003) and 
us also fitted a single elliptical Gaussian component to the innermost core for each epoch. 
The resulting major axis orientation is compared with the SPA orientation using the C1-C2 
position angle in Table 1. The differences are less than the uncertainties in the position 
angle determination, indicating that although the angular rcsohition of the 43 GHz maps is 
inadequate to determine details of the inner core substructure, the orientation of the inner 
core is probably accurately determined by Gaussian component modeling. 
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Table 1. Gaussian Model Parameters 



Epoch 


2 Component Gaussian 


EUip. Gaussian 


ASPA 


yr 


SPA(°) 


r(mas) 


SPA(°) 


AH 


1998.76 


207.2 ±9.3 


0.07 ±0.014 


200.0 ± 10.6 


7.2 ± 14.1 


1998.97 


191.5 ±8.3 


0.08 ±0.009 


190.1 ±7.1 


1.4 ± 10.9 


1999.16 


184.3 ±6.8 


0.08 ±0.007 


180.8 ±7.4 


3.5 ± 10.0 


1999.41 


193.0 ±8.5 


0.07 ±0.009 


194.3 ±9.2 


-1.3 ± 12.5 


1999.73 


180.3 ±9.1 


0.06 ±0.007 


170.9 ±7.9 


9.3 ± 12.1 


2000.01 


176.7 ±8.9 


0.09 ±0.007 


179.0 ±8.5 


-2.3 ± 12.2 


2000.31 


194.3 ±4.5 


0.12 ±0.007 


194.7 ±4.3 


-0.4 ±6.2 


2001.60 


188.9 ±5.3 


0.14 ±0.008 


185.2 ±5.2 


3.7 ±7.5 


2002.05 


186.9 ±6.8 


0.15 ±0.008 


187.5 ±6.9 


-0.6 ±9.7 


2003.82 


172.3 ±5.8 


0.09 ± 0.008 


175.0 ± 8.4 


-2.7 ±10.2 
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2.1. Model Component Uncertainty Estimates 

We used the program DIFWRAP (Lovell 2000) to estimate uncertainties in fitted Gaus- 
sian components. DIFWRAP is a 'front-end' graphical user interface to DIFMAP which 
allows the user to constrain a selected subset of model parameters over a fixed grid of values 
while allowing other parameters to freely vary. For each set of model parameters, the visi- 
bilities are self-calibrated and hybrid mapped, with the adjustable model components varied 
so ensure a best-fit to the visibihty data. This method has the advantage that it allows for 
the inter-dependence of model components to be taken into account during model fitting. 
The resulting self-calibrated maps are displayed for visual inspection with the value of 
each set of model components. 

We used DIFWRAP to vary the separation and SPA of the core subcomponents CI and 
C2 in steps of 0.02 mas and 5° respectively, followed by self-calibration and model fitting 
for each set of C1-C2 positions. We generated a sequence of hybrid maps corresponding 
to gridded values of separation and SPA over the ranges 0.02 to 0.20 mas and 180° to 240° 
respectively. At each epoch, the map corresponding to the minimum value was differenced 
from each of the other maps. In order to determine the uncertainties in the separation and 
SPA of the C1-C2 component pair, we adopted the heuristic criterion that the RMS noise 
level of the differenced map associated with a particular core separation and SPA did not 
increase by more than 50% from the lowest value for any map. 

Comparison of the SPA uncertainties of Stirling et al. (2003) with the present data 
shows that the former are significantly smaller than ours: The mean uncertainties are ±5.0° 
and ±7.3° respectively. This is surprising, since although the detailed method to determine 
uncertainties differs, the instrumental technique, mapping procedures, and resulting maps 
are nearly identical (e.g. Figure 1). More important, the published Stirling et al. SPA data 
does not agree with their proposed periodic model: The formal fit of SPA angle vs. epoch to 
the periodic function described in Stirling et al. (2003) can be rejected at 99.9% confidence. 
Hence, in order to compare our SPA data with those of Stirling et al. (2003) in a consistent 
manner, we solved for a constant multiplicative factor to apply to their SPA uncertainties 
which would result in a reduced = 1 when their model was applied. The resulting factor 
(1.5) made their uncertainties comparable with ours. In the following, we refer to the Stirling 
et al. dataset with larger uncertainties as the rescaled Stirhng et al. data set. 
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3. 



Results 



The variation of inner core SPA's versus epoch were compared with the periodic model 
of StirUng et al. as well as least-squares fitted periodic and constant models. We analyzed 
measurements derived from the ten epochs reported in this paper, the data of Stirling et 
al., and a combined dataset consisting of all 33 epochs. The fitted parameters as well as 
resulting reduced ^ and agreement factors are listed in Table 2. The periodic models have 
the form 



Column 1 lists the model being tested, while columns 2-5 list the model parameters. Columns 
6 and 7 list the SPA dataset and number of epochs. Column 8 lists the reduced of the 
model fitted to the data, while column 9 lists the corresponding goodness of fit probability 
(i.e., the probability that an experimental data set consistent with the assumed model would 
result in a reduced ^ value at least as large as the measured value) . We refer to this as the 
acceptance probability in the following. 

Figure 3a is a plot of core SPA (determined from C1-C2 position angle) versus epoch 
for all ten epochs fisted in Table 1, along with the null hypothesis model (non- varying SPA, 
solid line). The reduced of the fit is 2.76 (acceptance probability -p — 0.01). Hence, at 
99% (3(t) confidence the present SPA data indicates significant variations in the orientation 
of the inner core. Similarly, the one component elliptical Gaussian SPA data has a reduced 
X^ of 2.24 [p = 0.04), implying significant variability at 2a confidence. The Stirling et al. 
SPA data also indicate variability at high confidence {p < 10~^). 

Figure 36 shows the same data, but with a least-squares fit to a periodic model (solid line, 
period 13.1 yrs), and the model of Stirling et al. (dotted line). While neither model can be 
excluded, the fits are not particularly good: The acceptance probabilities are p = 0.14,0.10 
respectively. By comparison, the one component elliptical Gaussian SPA data has a best-fit 
solution with a 1.43 yr period and is in good agreement {p = 0.38). The large diff'erence in 
least-squares periodic solutions between the one and two-component SPA data, in spite of 
their similar values (Table 1), indicates that ten sampling epochs is probably insufficient to 
determine the true periodicity of the SPA data if it exists. 

Figures 3(c) shows combined SPA measurements from Stirling et al. (open circles) along 
with the present data (solid squares), and the Stirling et al. periodic model (dashed line). 
The reduced x^ of the fit is 1.88, or an acceptance probability p — 0.003. However, we have 
argued (see §2.1) that the uncertainties reported by Stirling et al. may be underestimated. 
Figure 3{d) shows the same data, but with rescaled Stirling et al. uncertainties. The Stirling 
et al. model (dashed line) and the least-squares best-fit periodic model (solid line) are also 




(1) 



-9- 



shown. The corresponding acceptance probabilities are 0.40 and 0.49 respectively. The best- 
fit period is 2.24 ±0.20 yrs, in excellent agreement with the 2.29 ±0.35 yr period determined 
by Stirhng et al. (2003). 
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Fig. 2. (a) 43 GHz map of BL Lac at epoch 1999.73. Contour levels increase in factors of 2 from 3.77 mjy per beam with 

a peak of 1510 mJy per beam. Dashed contours are negative. (V) Same as (a), but expanded to show only the inner core and 
centroids of best-fit Gaussian subcomponents CI, C2, and C3. The structural position (SPA) used to determine the orientation 
of the inner core is defined as the angle between subcomponents CI and C2. For this epoch, the SPA was 180.3° it 11.7° 
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Table 2. Model fits to BL Lac core structural position angles versus epoch. 



iViOtiei 




Model Parameters 








Xr 


PrnK 




A 

^0 


P{yr) 


to 












Stirling et al. 


12.2 


2.29 


1997.9 


191° 


This paper C1-C2 


10 


1.79 


0.10 












This paper Gaussian 


10 


1.47 


0.18 












Stirling et al. 


23 


2.28 


0.001 












Stirling rescaled 


23 


1.01 


0.44 












Combined 


33 


1.88 


0.003 












Combined rescaled 


33 


1.05 


0.40 


Best-fit periodic 


12.1 


13.1 


1996.3 


179° 


This paper C1-C2 


10 


1.61 


0.14 




12.0 


1.43 


1997.2 


183° 


This paper Gaussian 


10 


1.07 


0.38 




12.2 


2.29 


1997.9 


191° 


Stirling et al. 


23 


2.27 


0.001 




12.2 


2.29 


1997.9 


191° 


Stirling rescaled 


23 


1.01 


0.45 




12.1 


2.24 


1997.9 


191° 


Combined 


33 


1.82 


0.004 




12.1 


2.22 


1997.9 


191° 


Combined rescaled 


33 


0.99 


0.49 


Non-varying 








187° 


This paper C1-C2 


10 


2.78 


0.01 










192° 


This paper Gaussian 


10 


2.24 


0.04 










192° 


Stirling et al. 


23 


7.12 


< 10-^ 










192° 


Stirling rescaled 


23 


3.16 


< 10-^ 










191° 


Combined 


33 


5.35 


< 10-^ 










192° 


Combined rescaled 


33 


4.98 


< 10-4 
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Fig. 3. (a) BL Lac inner core structural position angle (SPA) vs. epoch for present observations (filled squares) and 

best- fit null hypothesis (non- varying) model (solid line). (6) Same as (a), but for best-fit periodic model (solid line) and original 
Stirling et al. (2003) periodic model (dashed line), (c) Combined SPA data of this paper (filled squares) and those of Stirling et 
al. (open circles) and the periodic model of Stirling et al. (dashed line), (d) Same as (c), but Stirling et al. SPA uncertainties 
have been scaled by a constant factor xl.5 (see text). The periodic model of Stirling et al. (dashed line) and best-fit periodic 
function (solid line) are superposed. 
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4. Discussion 

Do the present observations support Stirling et al.'s claim of core orientation periodicity? 
The results are ambiguous: While the SPA's derived from the present data are consistent 
with the Stirling et al. model, they are also consistent with several other periodic models. 
This is not surprising given the relatively large uncertainty of each SPA measurement with 
respect to the half- amplitude of the putative periodic variations. 

Given this ambiguity, it is useful to inquire whether agreement with the Stirling et al. 
model is statistically significant. The acceptance probabihty of the combined dataset to the 
best- fit periodic model {p — 0.49) is only marginally better than that of the rescaled Stirhng 
et al. dataset alone (p = 0.45). However, even with the relatively large uncertainties of the 
present SPA data, would a different set of measurements been able to falsify the periodic 
model of Stirling et al.? In order to test this, we randomly reassigned the epochs of our 
ten SPA measurements and calculated the goodness of fit to the Stirling et al. model. We 
generated 10,000 shuffled datasets and calculated the resulting reduced ^-nd acceptance 
probabilities for all cases. More than 96% of these shuffled datasets had acceptance proba- 
bilities p < 0.05, i.e., they could be rejected at 95% {2a) confidence or higher. We conclude 
that the current SPA data does provide independent confirmation (at 2a confidence) of SPA 
variations similar to those reported by Stirling et al. during the time interval spanned by 
both sets of observations. 

A more critical question is whether the combined dataset provides convincing evidence 
for (or against) periodic SPA variations. Since the combined dataset spans less than two 
periods with adequate samphng intervals (1998.23 - 2002.05), it is difficult to to make strong 
case for periodicity as opposed to episodic variations. If the variations are regular and 
periodic, measurement of SPA orientations over a longer timescale provides a critical test 
of periodicity. We can use the three epochs observed since the last SPA measurement of 
Stirling et al. (epochs 2001.60, 2002.05, 2003.82) to test periodic models. A goodness 
of fit analysis of these three measurements using the periodic model of Stirling et al. and 
the best-fit model in Table 2 results in acceptance probabilities p = 0.08 i.e., rejection of 
the model at shghtly less than 2a confidence. While this result is not definitive, it calls into 
doubt the periodic nature of the SPA variations. 

If the inner core of BL Lac is precessing, one might expect periodic variations in the radio 
fiux, since the relativistically beamed jet emission will be strongly modulated by the line of 
sight angle. The Doppler-boosted flux will depend on the Lorentz factor of the emission 
region and its angle to the observer's line of sight. Denn, Mutel, and Marscher (2000), in 
analyzing 17 epochs of BL Lac radio jet maps, fit a helical jet model inclined at an a angle of 
9° to the observer's line of sight. The four observed superluminal component Lorentz factors 
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were in the range 2.8 < 7 < 7.0 {Hq — 70 km s ^ Mpc ^). Since the observed half- angle of 
precession is 12°, the deprojected angle is ±1.8°. The expected flux ratio would be 

'S'min 

where 5 is the Doppler factor, a is the source spectral index, and 6max,min = 9° ± 1.8° are 
the extremal values of the line of sight angle. At centimeter wavelengths BL Lac is nearly 
always a flat spectrum source (a ~ 0), so the expected flux ratio over one precession cycle 
ranges from 1.1 (7 = 2.8) to 2.5 (7 = 7.0). Stirling et al. (2003) used their kinematic model 
for BL Lac and obtained similar results. 

Kelly, Hughes, Aller, & Aller (2003) have recently reported detection of quasi-periodic 
variations in the radio flux of BL Lac. They analyzed 20 years of synoptic radio flux mea- 
surements at 4.8, 8.4, and 14.5 GHz from the University of Michigan Radio Astronomy 
Observatory (UMRAO) using a cross-wavelet transform algorithm. They report significant 
frequency-dependent peaks: 1.4 yr (4.8 GHz), 3.7 yr (8.4 GHz) and multiple periods (0.7, 
1.6, and 3.5 yr) at 14.5 GHz. These are near, but not definitely not within the uncertainty 
range of the proposed precession period of 2.29 ± 0.35 yr of Stirling et al., or the amended 
period 2.24 ± 0.20 yr determined using all 33 epochs. Villata et al. (2004) have searched for 
periodicities in historical radio and optical light curves for BL Lac from 1968-2003. They 
analyzed mult i- frequency radio observations from 4.6 to 37 GHz using several statistical 
methods suitable for unevenly spaced samples. Although they found that radio outbursts 
tend to repeat every ~ 8 years, there is no evidence for fiux modulation at the proposed 
precession period. 

The non-detection of flux periodicity at the precession period is not necessarily a serious 
objection to the precession hypothesis since a large fraction of the total radio flux arises from 
the extended jet which clearly does not systematically vary in position angle (Denn, Mutel, 
and Marscher 2000; Stirling et al. 2003). 

5. Summciry 

We have searched for conflrmation of the periodic inner-core orientation changes re- 
ported by Stirling et al. (2003) in the parsec-scalc radio jet of the AGN BL Lac. We an- 
alyzed 43 GHz VLBA observations over ten epochs from 1998.76 - 2003.82 which overlap 
the time range reported by Stirling et al. As expected, the resulting maps are very similar 
in quality and angular resolution to those of Stirling et al. Using a Gaussian model-fitting 
procedure similar to that of Stirling et al. (2003) but with different criteria for estimating 
the uncertainties, we find: 



5 (7> ^max) 
5 h,Omin) 
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1. We confirm (at 95% confidence) variabifity of the structural position angle (SPA) of 
the the inner core. The acceptance probability of the null hypothesis (no variation) 
was p — 0.04, 0.01 for the one and two-component Gaussian model fits to the inner 
core. 

2. Our SPA orientation measurements are consistent with those of Stirling et al. (2003) 
during overlapping epochs, after multiplying their uncertainties by a constant factor of 
1.5. This adjustment also makes their periodic model agree (in a statistical sense) 
with their SPA data: The acceptance probability is reduced from p < 10~^ iop — 0.44. 

3. Our SPA data is best- fit by periodic functions with periods 12.0 yr (p = 0.38) and 12.1 
yr (p = 0.14) for one and two-component Gaussian models respectively. However, the 
2.29 yr periodic model of Stirling et al. was also a plausible fit to our data (acceptance 
probabilities p = 0.18, 0.10). 

4. Combining the rescaled Stirhng et al. (2003) SPA data with the present dataset, we 
find a best-fit periodic model with a period of 2.24 ± 0.2 yrs (p = 0.49), very close to 
the 2.29 year period found by Stirling et al.. 

5. Using the best-fit model, a test of agreement with three SPA measurements made 
since the last Stirhng et al. (2003) epoch (2000.6- 2003.78) results in an acceptance 
probability p = 0.08, i.e. rejection of the periodic model at slightly less than 2(7 
confidence. 

If the relativistic jet is precessing, the consequent variation in line-of-sight angle implies 
periodic flux modulation of the inner core component. We estimate the fiux ratio to be 1.1 

to 2.5 based on jet parameters modeled from studies of superluminal components (Denn, 
Mutcl, and Marschcr 2000). However, studies of the radio and optical light curve of BL Lac 
over several decades (Kelly, Hughes, Aller, & AUer 2003, Villata et al. 2004) find no periodic 
variation at the suggested precession period. 

In summary, our although our SPA data provides confirming evidence for orientation 
variability in the inner core of BL Lac consistent with the measurements of Stirling et al. 
(2003), the evidence for periodic variations is less compelling. A definitive resolution will 
require continued monitoring of the inner core radio structure over several periods. 

We thank Casey Dreier for assistance with the data reduction. 
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